arXiv: 1507.04986vl [math.AP] 17Jul2015 


FRACTIONAL DISCRETE LAPLACIAN 
VERSUS DISCRETIZED FRACTIONAL LAPLACIAN 


OSCAR CIAURRI, LUZ RONCAL, PABLO RAUL STINGA, JOSE L. TORREA, AND JUAN LUIS VARONA 


Abstract. We define and study some properties of the fractional powers of the discrete Laplacian 

(~A h ) 3 , on 2^ = AZ, 

for h > 0 and 0 < s < 1. A comparison between our fractional discrete Laplacian and the discretized 
continuous fractional Laplacian as h —> 0 is carried out. We get estimates in £°° for the error of 
the approximation in terms of h under minimal regularity assumptions. Moreover, we provide a 
pointwise formula with an explicit kernel and deduce Holder estimates for (— A^) 3 . A study of the 
negative powers (or discrete fractional integral) (—A^) —s is also sketched. Our analysis is mainly 
performed in dimension one. Nevertheless, we show certain asymptotic estimates for the kernel in 
dimension two that can be extended to higher dimensions. Some examples are plotted to illustrate 
the comparison in both one and two dimensions. 


1. Introduction and main results 


The fractional Laplacian, understood as a positive power of the classical Laplacian, has been 
present for long time in several areas of Mathematics, like Fractional Calculus and Functional Analysis. 
However, although this operator was used for some differential equations in Physics, it was not until the 
first part of the last decade when it became a very popular object in the field of Partial Differential 
Equations. Indeed, equations involving fractional Laplacians have been one of the most studied 
research topics in the present century. Nowadays hundreds of papers can be found in the literature 
with the phrase “fractional Laplacian” in their titles. We recall that the fractional Laplacian on R ra 
can be defined for 0 < s < 1 and good enough functions u as, see ( 121 . 


(—A ) s u(x) = c n , s P.V. 


u(x) - u{y) 

\x — y\ n + 2s V 


The trigger that produced the outbreak in the field was the paper by L. Caffarelli and L. Silvestre 
[3j. The essence of that work is to immerse the fractional Laplacian (a nonlocal operator) on R n into 
a (local) partial differential equation problem in K" +1 , where the known methods of PDEs can be 
applied. Since the appearance of that paper there has been a substantial revision of a big amount 
of problems in differential equations where the Laplacian is replaced by the fractional Laplacian. Of 
course in this trend the Numerical Analysis should have to be present. Several recent works can be 
found in this direction, see [a [a nu nn Eg. One of the main troubles to overcome in the numerical 
approach is the nonlocality of the operator. 

On the other hand, in the past few years the language of semigroups has been used as a versatile tool 
in the study of different aspects and scenarios for the fractional Laplacian with a great success. See 
mii using m mg. In particular, in [5] it was possible to define the fractional discrete Laplacian 
applying this approach. Therefore the semigroup language gives us the chance to define and study 
the fractional discrete Laplacian. An obvious question arises immediately: can this language be used 
to develop some kind of Numerical Analysis for the fractional Laplacian? This is exactly our aim 
here, to present a discretization method to approximate the fractional Laplacian. 
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Let us give a brief description of our results. Consider a mesh of size h > 0 on R. given by 
Z/j = { hj : j £ Z}. The discrete Laplacian A h on the mesh is given by 

~A hUj = l - 2uj 

Here Uj, j £ Z, is a function on Z/j. The fractional discrete Laplacian (—Ah) s is defined with the 
semigroup language: 

1 poo u 

( L1 ) (- A h ) s u 3 = J o (e tAh Uj - Uj ) 

where Wj{t) = e tAh Uj is the solution to the semidiscrete heat equation 

j d t Wj = A h Wj, in Z h x (0, oo), 

[ Wj(0) = Uj, on Z/j. 

Now given a function u = u(x) defined on R, we consider its restriction r^u to the mesh as defined by 
( rh‘u)j = rh‘u(j) := u(hj). Our goal is to estimate, as a function of the size h of the mesh, differences 
of the type 

|| (—A h ) s (r h u) - rh({—A) s u) || <00 . 

Certainly the estimates will depend on the smoothness of the function, as our first main result shows. 

Theorem 1.1. Let 0 < /? < 1 and 0 < s < 1. 

(i) Let u £ C' 0 ’^(R) and 2s < (d. Then 

|| (-A h) s (rhu) ~ rh((~A) s u) \\(°° < C[u] G o,e(K)h^~ 2s . 

(ii) Let u £ C' 1,/3 (R) and 2s < fd. Then 

\\D + (-A h ) s (r h u) - r h (£(-Ayu)\\poo < C[u\ c i,p (V) h^ 2s . 

(in) Let u £ C' 1 ’ /3 (R) and 2s > fd. Then 

||(—A h ) s {r h u) - r-^((-A) s u)||^oo < C'[m] c -i./3 (r) /i / 3_2s+1 . 

(iv) Let u £ C k,/3 (M.) and assume that k + fd — 2s is not an integer, with 2s < k + (d. Then 
\\D l + (-A h y(r h u) -r h (| r (-A) i u)l|/~ < C[u] c ^ m h^ 2s+k ~ l , 

where l is the integer part of k + fd — 2s. 

The constants C above depend only on s and fd, but not on h or u and j- denotes the standard 
derivative. For a function Uj on the mesh, we let 


D + Uj = 


Uj+l Uj 


and by D l + we mean that we apply the operator D + to Uj l-times. 

To prove this result we need the explicit formula for (—A/,) s presented in the following statement. 


Theorem 1.2. Let 0 < s < 1 and £ s := 
(a) For 0 < s < 1 and u £ £ s we have 


: Z;, 


II || \ s I a m | 

— 2^ (i + | m |)i+2« < 00 
mez v 117 


( 1 . 2 ) 


(-A h ) s Uj = ( u o - u m)K s {j - m), 


m^j 


(1.3) 


where the discrete kernel K s is given by 

4 s r(l/2 + s)T(|m| — s) 


KJm) = 


V^r|r(-s)|r(|m| + 1 + s) ’ 


m £ Z \ {0}. 


(b) If u £ £q then lim (—Ah) s Uj = Uj. If u £ £°° then lim (—A h) s Uj = —A h'Uj. 

s— > 0 + s —►!- 
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(c) For 0 < s < 1 there exists a positive constant C s such that 


(1.4) 


K s (m) < 


C a 

|m| 1+2s ' 


Observe that, in view of (1.2), the fractional discrete Laplacian is a nonlocal operator on Z^. Notice 
also that our definition is neither a direct discretization of the pointwise formula for the fractional 
Laplacian, nor a “discrete analogue”, but the s-fractional power of the discrete Laplacian. 


Remark 1.3 (On the negative powers). We recall here the nice work by E. M. Stein and S. Wainger, 
see |22| , in which they consider a fractional integral type operator on Z given by convolution with the 
kernel |?ti| - ( 1-2s ). It is interesting to notice that if we define, again with the help of the semigroup 
language, the negative powers (or fractional integrals) (—A^) _s then the kernel K_ s of this operator 
satisfies K_ s {m) < C s |m|~ < ' 1 ~ 2s ). Moreover in the case 0 < s < 1/2, there exist positive constants c s 
and C s such that |A'_ s (m) — c s |to|“( 1_2s )| < C s \m\~^ 2 ~ 2s \ see Sectionji] 

Remark 1.4 (Probabilistic interpretation). Let u be a discrete harmonic function on Z^, that is, 
—A hU = 0. This is equivalent to 

1 1 

u i ~ 2 U 7+ 1 + 2 Uj ~ 1 ' 

This mean value identity says that a discrete harmonic function describes the movement of a particle 
that jumps either to the adjacent left point or to the adjacent right point with probability 1/2. 
Suppose now that Uj is a fractional discrete harmonic function, that is, (—A h) s Uj = 0. Then we have 
the following mean value identity: 

Uj = =- ^2 u m K s (j - to), 


where £ s = o ^ s( m ) = 2 ^p 1 ^ 2 / 1 /^ ■ In a parallel way we understand this last identity as saying 
that a fractional discrete harmonic function describes a particle that is allowed to jump to any point 
on Z/j (not only to the adjacent ones) and that the probability to jump from the point j to the point 
to is / K s (j — m). As s —> l~ the probability to jump from j to a non adjacent point tends to zero, 
while the probability to jump to an adjacent point tends to one, recovering in this way the previous 
situation. As s —> 0 + , the probability to jump to any point tends to zero, so there are no jumps. 


Formula (1.2) makes it possible to prove Holder estimates parallel to the corresponding estimates 
proved for the classical Laplacian in (2U1 . The result we get is the following. For the definition of 
discrete Holder spaces C^’ /3 see Section Klj 


Theorem 1.5. Let 0 < /3 < 1 and 0 < s < 1. 

(i) Let u € C° ,/3 and 2s < /3. Then (—A h) s u £ C®’ 13 " s and 


II ( — Aft) 11 £.0,3-2 b < Cj|u|| c o,3. 

(ii) Let u £ C^ 13 and 2s < fd. Then (—A h) s u £ C^ ,/3-2s and 

||(-A^) S u|| c i, 3-2, < C||u|| c i,3- 

(in) Let u £ C^’ 13 and 2s > fd. Then (—A h) s u £ (j^P~ 2s+ i an d 

||( — A/ 1 ) S U|| C ,0,3-2B+1 < C\\u\\ c y. 


(iv) 


Let u £ C^’ 13 and assume that k+/d—2s is not an integer, with 2s < fc+/ 3. 
where l is the integer part of k + fd — 2s and r ) = k + fd — 2s — l. 


Then (- A h ) s u £ C l f 7 , 


The constants C > 0 above are independent of h and u. 
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Although the proof of Theorem 1.1 i| is not trivial, one could say in a naive way that such an 
approximation result is in some sense announced by the results in Theorem |1.5[ Indeed, the fractional 
discrete Laplacian maps into C^~ 2s . The continuous version of this property is also true for the 
fractional Laplacian, so the restriction of (—A ) s u to the mesh Z^ is in C^~ 2s whenever u £ C a (M). 

So far the results shown above were developed in the one dimensional case. Analogous results to 
Theorem 11 . 2| and Theorem |1.5| can be extended to higher dimensions. These multidimensional results 
will appear elsewhere. Our techniques depend strongly on the explicit expression for the kernel of 
the fractional operators. Nevertheless, we cannot find such explicit expressions for the kernels in the 
two (or higher) dimensional case. This yields technical difficulties that we cannot overcome to prove 
a multidimensional version of Theorem EH However, we have been able to get asymptotics for such 
kernels in dimension two that we believe could also be reproduced for the multidimensional case. 

At the end of the paper we will show some pictures of examples drawn with Mathematica that 
illustrate our fractional discrete operators, both in one and two dimensions. We take examples of 
functions u in the continuous variable for which / := (—A ) s u is explicitly known. We plot / together 
with the discrete function fj = (—A h) s (,i"hu)j. We also consider some explicit solutions u = (—A )~ s f 
for given /. We compare the graph of u with the solution Uj = (—A h)~ s { r hf)j- In this regard, we 
are addressing the question of whether the Poisson problem for the fractional Laplacian 

(-A ) s u = f, inr, 

can be discretized by using our formulas, and whether the solutions to this discretized problem 
converge in some sense to the solutions of the continuous Poisson problem. 

Several authors have been interested in solving the fractional Poisson equation in a numerical or 
discrete way. We mention here the recent works rnrnmm, see also the references therein, for the 
numerical approach, and l29l for a one-dimensional discrete approach. 

As we remarked earlier, the strategy used to obtain our results is the language of semigroups. Since 
the semidiscrete heat semigroup is given in terms of modified Bessel functions, see Section [2] we will 
use exhaustively some properties and facts about these functions. 


The structure of the paper is as follows. Sectionals devoted to the proofs of Theorems 1.2 and 1.5 


The proof of Theorem 0 is included in Section [3] In Section [4] we analyze the negative powers 
(—A h)~ s - The asymptotics for the kernel of the fractional discrete Laplacian and fractional integral 
in dimension two are studied in Section [5j In Sections [6] and [7] we show our pictures. Finally, some 
technical proofs needed on the way and the properties of the Bessel functions are collected in Section[8] 
By C s we mean a positive constant depending on s that may change in each occurrence, while by C 
we will denote a constant independent of the significant variables. 


2. The fractional discrete Laplacian: proofs of Theorem 11.21 and Theorem 11.51 

Given Uj : Zj, -> K, the solution to the semidiscrete heat equation with initial discrete temperature 
Uj can be written as 

(2.1) e tAh Uj = ^2 G(j - m, ^s)u m = ^ G(to, f€ [ 0,oo), 

raEZ mEZ 

where the semidiscrete heat kernel G is defined as 

G(m, t) = e~ 2t Im{2t), m £ Z. 

This follows by scaling from the case ft. = 1 of isms. Here, /„ is the modified Bessel function of order 
v whose properties are collected in Section [8] By (8.5) and (8.6) the kernel G(m,t ) is symmetric in 
m, that is, G(m,t) = G(—m,t), and positive. 

Proof of Theorem\l.£\ First we will check that e tAh Uj is well defined. Indeed, if N > 0, by using the 


asymptotic of the Bessel function for large order ( |8.11| ), then 

„ — 2 1 


H G(m, ^)\uj- m \ < Ce 2t 

|ra|>JV \m—j\>N 


f m \ (1 -)- 1 777, — j|) 1 ' t ' 2s 


r(|m| + 1) (l+|m-j 


U Vki^<o tM \\n\u s . 












(a). Define 

( 2 . 2 ) 
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K s (to) = 


G(m, t) 


dt 

fi+s 


\n-s)\ Jo 

for m/0, and K s { 0) = 0. The symmetry of this kernel in m follows from the symmetry of G(m,t). 
Therefore it is enough to assume that m £ N. To get ( |1.3[ ), we use (8.12) with c = 2 and v = m. On 
the other hand, it is easy to show (see for example [6] for the case h = 1) that e tAh l = 1. Hence, 
from and ( |2.l| ), 


i r°° _ 

(-A h y Uj = ——l 53 G(j-m,^)(n m -w,) — 

V ' 0 m^j 

1 /*°° 

= f(Z^) 53 “ u j) J o G 0 ~ m ’ w) 


m^j 


dt 

t 1+s 


1 

v—d 53 I G (J ~ m ’ r ); 

' '' m^j ® 


/i 2s |r(- 

53 ~ u m)K s (j - m). 


dr 

>1+S 


h 2 




For the interchange of summation and integration in the second equality, we consider the two terms 

„ i. t \ i | dt , , f°° \ ^ t \ dt 

m^j 


53 G l? - 

m^j 


By using (1.4) we see that the first term above is bounded by C s YJm+j \ m ~ J 
finite for each j because u £ t s . For the second term we use again ( |1.4[ ). 

(b). Suppose first that u £ £q. We have 


t i+s- 

— '!|-( 1 + 2s ) \u m \, which is 


h 2s (-A h ) s Uj = Uj 55 K s{j ~ m) - 53 K s{j ~ m)u m =: UjT ± - T 2 . 

m^j m^j 

We write 7) = T ltl + T\ 2 , where 


T b 2 =rrCT53/ rrraE / G ( TO ^ 


dt 


1 


l r (-s)l 




* 1+s |r(-s)| 


ra^O 1 


dt 

H+s' 


We are going to prove that Tip and T 2 tend to zero, while T ij2 tends to 1, as s —> 0 + . Let us begin 
with T 12 . By adding and subtracting the term m = 0 in the sum and using (8.7), we get 


t i,2 = 


l r (-«)l 


*/o(2t) 


t l + S 


dt . 


By noticing that |T(—s)|s = F(1 — s) and that, by (8.10), we have 


1 

|r(-s)| 


2t I 0 2t , G r°° 

; dt < T——TV / t 

t i+s - |r(-s)| A 


- 1 / 2 - 1 - dt = 


c 

l r (-s)|(l/2 + s) : 


we get T 1i2 —> 1 as s —> 0 + , as desired. Next we handle the other two terms Tip and T 2 . On one 

1 


hand, by (8.9), 


T 


1,1 


— y 

(-A 


1 


-2 1 


l r (-s)l r(|w| + “'O 


fl m l 


1 


- iri-sd 53 


i 


dt 

t 1+s ~ |r(-s)| ^ o r(|m| + l) |m| —s’ 


and the last quantity tends to 0 as s —► 0 + . On the other hand, for T 2 , we use (1.3) to obtain 

\rr i ^ G « r(|j — m\ - s) 

2 “ \n-s)\^n\j- m \ + i + S ) 1 


The constant G s remains bounded as s —>• 0 + . Since u £ £q, by dominated convergence, we have that 
the sum above is bounded by ||n||^ 0 , for each j £ Z. Therefore T 2 tends to 0 as s -A 0 + . 
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For the second limit, suppose that u £ . By using the symmetry of the kernel K s we can write 


where 


h 2s (-A h y Uj = s 1 + s 2 , 

Si = K s (l)(-u j+ i + 2 Uj - Uj- 1 ), and S 2 = ^ K s (m)(uj - uj- m ). 

|m|>l 


Next we show that K s (l) —> 1, while S 2 —> 0, as .s —> 1 , which would give the conclusion. By (1.3) 
we have 

lim K,( 1) = lim 4 ^, 1 7 , \ r(1 / 2 + l > = Urn + = !. 

s-n- y l- v^r|r(-s)|r(2 + s) 8 ->i- + s) ^(3) 

On the other hand, by ( |1.3| ), S 2 is bounded by 

„ 4 s r(i+s) ^ r(lml-s) 
ll U ll^ $5||u||/o. / E r(H + l + s)’ 

|m|>l 1 V n |m|>1 vl 1 ' 

which goes to zero as s —► 1“. 

(c). Observe that the estimate in (1.4) follows from (1.3) and Lemma 8.2 


□ 


Consider the first order difference operators on Z/, given as (observe that we already introduced 
some definitions in the statement of Theorem 0 >- 

D +Uj := Uj+ \~ Uj , and D_u 3 := Uj . 

h h 

For 7,77 € No, we define Dy v _Uj := DlDfuj. Here, by D± we mean that we apply k times the 
operator D± and D*± is the identity operator. 

Definition 2.1 (Discrete Holder spaces). Let 0 < /3 < 1 and k £ No- A bounded function u : Z^ — > ffi. 
belongs to the discrete Holder space C 1 ^' 13 if 


[D+'Lu] c o,i3 := sup 


| D™Uj - D™i 


< C < oo 


h,P\j — m\P 

for each pair 7 , 77 £ No such that 7 + 77 = k. The norm in the spaces is defined in the usual way: 
||7t|| r fc,/3 := max sup \D'y r, _u m \ + max [D'y r, _u] r o,f>. 

h J+V<k me z h ’ 7 +V=k ’ h 

Proof of Theorem m By iteration of (i), (ii) and (Hi) we can get (iv). Let us start with the proof 
of (i). Let k,j £ Z. We have 

(2.3) k-a h ) s u k - (-A h y Uj \ = 1 L\s 1 + s 2 \ 


where 

(2.4) 


h 2s 


Si := ^2 ( U k - Uk+m - Uj + Uj + m)K s (m ), 


and S 2 is the rest of the sum over \m\ > \k — j |. By (1.4) 

\m' R 


Si<C s [u\ c o,/>h 0 


|l+2s 


< C a [u] c o.fih f> \k-jf- 2a . 


l<|m|<|fc-j| 


For S 2 we use that \uk — Uj\ < [v\ c o,ph^\k — j\& and (1.4) again to get 

S 2 < C s [u] c eh 0 \k - jf J2 HI" 1 " 2 ' < C a [u] c o,„h f, \k - jf 2 ~ 2a . 

\m\>\k-j\ 

For (ii), it suffices to prove that D±(—Ai l ) s u £ C° ,/3 ~ 2s . For this it is enough to observe that D± 
commutes with (—A^) s and then use (i). 
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For (in) we are going to use (2.3). Without loss of generality, take m £ N. We split the sum in 


(2.3) by taking the terms u k — u k +m and Uj — Uj+ m separately. The following computation works for 
both terms, so we do it only for the first one. It is verified that 


(2.5) 

Therefore, 

( 2 . 6 ) 


m— 1 

Uk-\-m Uk = h ^ ^ 7 * 

7=0 


m— 1 


Uk 


Uk+m = [hmD + u k - h E D + u k+ ^j - hmD + u k ■ 


7=0 


On one hand, by taking into account that the kernel K s (m) is even, we get 

( 2 -7) E ( hmD + u k )K s (m ) = hD + u k El mK s (m ) = 0. 

On the other hand, since u £ CE> the first term in the right hand side of (2.6) is bounded by 

771—1 


771—1 


(2.8) h E \D+u k - D + u k+1 1 < h 1+l3 [u] c i,f< £ M i h^luU., |m|^|m| = [u] c i,ft(h\m\) 1+l3 . 

7=0 7=0 


Pasting (2.7) and (2.8) (and their analogous for Uj — Uj +m ) into (2.41, we conclude that 


IE < C s [u] c i,f>h 1+ P E 


\m 


1+/3 


I l+2s 


< C s [u\ c i,f,{h\k - j\) 1+f) - 2s . 


Now we deal with S' 2 . By ( |2.5[ ), 

| (^fc Llj ) ( a fc+m ^j+m)! | (^(fe— j)+j ) (^(fc—^i+m)| 


k—j-1 


— h ^ 1 — (U ] & A h> ^ 17Ti | ^ | /c j |. 


7=0 


Hence, 


I-S 2 I < C[u] c i,ph 1+ P\k — j\ E \m\ p K s (m) < C[u\ c i,ph} + P\k - j| 1+/3 


□ 


3. Error of approximation: Proof of Theorem 11.11 

Recall that a continuous real function u belongs to the Holder space C^ER) if u £ C fe (R) and 

(k) \u( k \x) - (y)\ 

u (fe) C'O^(R) := sup J -;- r-o - - < OO, 

z.ySR \X-y\P 
x^y 

where we are using u^ to denote the /c-th derivative of u. The norm in the spaces C fe,/3 (R) is 

k 

IMIo^R) := E ll u EU°°(R) + [ u ^]c°'0(R)- 

l -0 


In order to prove Theorem 1.1 we need a preliminary lemma. Let us define the following constant: 

4 s r(l/2 + s) 


(3 - 1} Cs 7r 1 / 2 |r ( -s)| 

Lemma 3.1. Let 0 < s < 1. Given j £ Z, we have 

(3.2) - ^ dy Ks{m) 


> 0 . 


'I y 


-h(j+m)\<h /2 \hj - y\ 1+2s h 2s 


< 


C, 


h 2s m 2+2s 


, for all m £ Z, 
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(3.3) 

and 

(3.4) 


dy 


< 


C s 


l\y-h(j+m)\<h /2 \hj-y\ 1+2s h 2s m 1+2 - 

hj - y 


for all m £ Z, 


E 


^\y— h U+ m )\ < h/ 2 \hj y l 1+ "‘ 

Proof. The change of variable hj — y = hz produces 
dz K s ( m ) 


■dy = 0. 


h 2s J\z-m\<i/2 M 1+2s 


h 2s 


< 


C s 

h 2s 


By using the Mean Value Theorem, 

1 


' \z—m |<l/2 


'\z— m|<l/2 

l 


l 


l 


\z\ 1+2s \m\ 1+2i 


dz 


,-2s 


TO 


l+2s 


- K s (to) 


UH+2S | m |l+2s 


dz 


<C S 


dz 


while by Lemma 8.2 


|l+2s 


- K s (m) 


' \z—m\ <1/2 l TO 

Cs 


2+2s 


c fl 


Ira 


2+2s ’ 


| to | 2 +2s ' 


Thus (3.2) follows. For (3.3), it is easy to see that 

1 dy <cj 


dy 


C s 


\y-( h ( j+m ))\<h/2 I hj - y\ 1+2s J\y-(hU+m))\<h/2 IM 1+2s h 2s \m\ 1+2a ' 

dy = 0. 


Finally, let us prove (3.4). By symmetry, we have 

/' {hj - y ) 


J\y-hj\<h /2 \hj-y\ 1+2 ' 

Moreover, by changing variables hj — y = z, we get 

z 7 v-^ f z 


E 


1 ,^ Z , '\ z ~ hrn \ <h / 2 I z 

m^O 


1+2 s 


d * = E 


* z-\-hd\<Zh j *1 

£^0 


l+2s 


d 2 = E 


J\u—h£\<h/2 


—U 

l+2s 


du, 


and the conclusion readily follows. 

Proof of Theorem o (i). We write, for each j £ Z, 

r u(hj) - u(y ) 

m& zJ\y- h U+™)\ <h /2 \hj ~ y\ 1+2 ‘ 

u{h{j + to)) - u(y) 
~z J \y- h (j+ m )\ <h / 2 \hj — y\ 1+2s 


r h ((-A) s u). = c a Y 

iez 

E 


dy 


= Cs 


dy 


+ Y ( u (hj) - u{h{j + to))) f 


dy 


m£Z 

mz/z 0 


=:c s {S 1 + S 2 ). 

By using that u £ C' 0,/3 (K) and (3.3), we have 


I Si I < C'[u] c o,p( M ) Y. 


\y—h(j+m)\<h/2 I hj y\ ~*~ 


h d dy 


ra£ Z 1 


\y— h{j+m)\<h/2 I ^3 V I 


< C s [u]co.n(s.)hd ’Y, 


1 




/i 2 s |m| 1+2s 


— C' s [ti]c’0,^( R )/l /3 


□ 
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Now we compare S 2 with (—A h) s (i"hu)j. Since u € C' 0, ^(K), by Lemma 


3.1 


we can see that 


Y ( u ( h 3) - HHj + m))) / 

mfZ J \y 


dy 


mG 

m^O 


-h(j+m)\<h /2 \hj~y\ 1+2s 


< Y \u(hj) - u{h{j + m))\ c s / 

-'ll 


- (-A h ) a {r h u)j 


dy 


K s (to) 


mGZ 

ra^O 


< C' s [u]c°.£(R) ^ 


mGZ 

m /0 


/i 2 s |to| 2 + 2s 


|y — ii(j+m)|</i/2 dj y | h 


< C s [u]co,n(mh& 


(ii). Observe that j- commutes with (—A) s and D + with (—A/,) s . Then 
||D+(-A ?l ) s (r /l u) -r^(£(-A) s u)|| < , oo 

< ||(-A/,) 8 £>+(r h u) - (-A^) s (r h ^ M )|| foo + ||(-A h ) 8 (r h ^ti) - r h ( A(-A) s m) || /bo . 
For the second term we just apply (i). As for the first one, by using the Mean Value Theorem, 

| (-A h ) s D + (r h u)j-{-A h ) s (r h £u) | 


h 2i 


Y K s (m) 


m^O 


w(/i(j + 1)) - u(hj) 


- u'(hj) 


u(h(j + m + 1)) - u(h(j + mn))) 


— u (h(j + to)) 


^ Y K s(m) 


ra^O 


(«'(&) - u'(/lj)) - (w'fe+m) - u'(h(j + to))) 


< C, 2s ]T K s (m)h p < Ch? 


— 2s 


m^O 


where ^ is an intermediate point between hj and h(j-\-1 ) and analogously In the last inequality 

we have used the hypothesis on the regularity of u. 


(in). By taking into account (3.4), we can write 

u(hj) - u(y) - u'(hj)(hj - y) 


r h ((-A) s u) = c s Y \ 

\y-h(j+m)\<h/2 


dy 


I hj - y | 1+2s 
u(h(j + to)) - u(y) - u'(hj)(h(j + to) - y) 


^^j\y~ h (j+m)\<h/2 \ h 3 V\ + 

+ Y (H h i ) ~ HHj + m )) _ u '( h J)( h J - Hj + m ))) / 

«/k 


dy 


dy 


mGZ 

m^O 


— ^.(j+m)|</i/2 2/1 


=:c s (S 1+ S 2 ). 


By the hypotheses and (3.3), we have 


ISil < CsMcijjb) Y \hmfh 


1 


ra£Z 


h 2s \m\ 1+2s 


— C s [m]ci,/3( R )/i 1+/3 2s . 


We compare c s S 2 with (—A/ l ) s (r/ 1 rt)j. Since K s (m) is even in to, we can write 

(-A h) a {r h u)j = Y (HH) - u{h(j + m)) - u'{hj)(hj - h{j + to))) • 


mGZ 

m /0 


Then (3.2) and the hypothesis on u give the result. 
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(iv). The proof in this case follows as in (ii) by iteration l times. 


□ 


4. The negative powers or fractional discrete integral 

Analogously to the fractional discrete Laplacian, the negative powers of the discrete Laplacian, 
also called the fractional discrete integral , is defined for a discrete function / : Z^ —> R. and jgZ as 

dt 


(-A h )~ s f 3 = 


1 

W) 


t A. 


u 


3 t l ~ s ' 


s > 0. 


By writing down the semidiscrete heat kernel and using Fubini’s Theorem we readily see that this is 
a convolution operator on Z: 

(-A h )~ s fj = h 2s E K_ s (j - m)f m , 
m {EZ 

where the discrete kernel I\- s is given by 

A '- (ra) = F hi G{m - t) J L -‘ 


It is worth to compare this formula for 7\_ s with the one for K s , see (2.2). 

Theorem 4.1. Let 0 < s < 1/2. We have the explicit expression 

(4.1) + K-.(O) = "“> 


\fit r(s)r(|m| +1 - s) 

Moreover, there exists a positive constant C s such that, for m £ Z \ {0}, 


- s) 


(4.2) 


K- S {m) - 


11 —2s 


< 


C' 


Proof. One can compute the integral defining the kernel explicitly, for 0 < s < 1/2, just by using 
formula (8.12). The estimate in (4.2) follows from (4.1) and Lemma 8.2 □ 


As we mentioned in the Introduction, a discrete analogue of the usual continuous fractional inte¬ 
gration can be defined for / : Z” —)■ K and 0 < A < n by 

1 


h /, = E 


mGZ n 

ra^O 


\m-J\ 


-Txfm, 


E. M. Stein and S. Wainger obtained £ p —> £ q estimates for this operator in [52]. In view of our 
expression for the kernel AT_ S in ( |4.1[ ) and (8.2 1 of Lemma 8.2 we notice that the one-dimensional 
version of the operator /i_ 2 S controls, up to a factor of h, our fractional discrete integral (— Ah)~ s . 
Therefore, mapping properties for (— Ah)~ s can be deduced from the results in [221 . 

5. Asymptotics for the two dimensional fractional discrete Laplacian 

In this section, we present an expression for the two dimensional kernel K s (m), m = (mi, m-f) € Z 2 , 
when both |mj| —¥ oo. We denote 11m.11 2 := (|mi | 2 + Im^l 2 ) 1 / 2 , for m £ Z 2 . 


Lemma 5.1. Let 0 < s < 1. Let K s be the kernel in (2.2) in dimension two. Then 
(5.1) 0 < K s (m) = ——I - hiqher order terms, when both —> ± 00 , * = 1 , 2 , 

v / sv ; ||to|| 2+2s 

where 


C2,s = 


4 s T(l + s) 

7r|T(-s)| 
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Proof. The ideas in this subsection are inspired in some results and techniques from [5]. First we list 
the ingredients that will be used. 

In 1151 Chapter 4, (5.05)] we find the following expansion for the ratio of two Gamma functions, 
for 2 £ R, z —> oo, 


1 + 


(a — b)(a + b— 1 ) 1 fa — b 


2z 


12 


^ (3(a + b ) 2 — 7a — 56 + 2) —— ^ + E[a, b, z ), 


(5. 2 ) = *«-• 

r(2 + b ) 

where E(a , 6 , z) is an integral (convergent if 3 + b — a > 0) whose absolute value is bounded by 
z ~(3 +b-a) times a constant depending on a and b. 

We will use formula (8.12 1 with c = 2, a = —s (where 0 < s < 1), and v = where rrii £ N, 
m, ^ 0. Indeed, we have 

( 5 . 3 ) f < 20 .-'* = 4 * r (l/2 + »)r(mi - ») 


/*oo 

/ e~ 2t I mi {2f)t 
Jo 


Let 

(5.4) 

Consider the Mellin transform 

(5.5) 


r/n F (m, + 1 + s) 

frrii if) = e~ 2t I mi (2t). 

dt 


/m, (s) = / frrii (0 


fl+s 


which is absolutely convergent, by (5.3), for 0 < s < 1. By the Mellin inversion theorem we can 
recover f mi ( t ) via the inverse Mellin transform, see for instance [26], that gives 


(5.6) 


pc+zoo 


2ni 


f mi {s)t s ds , 


which is well defined for c a real value between 0 and 1 ; the integral above is a line integral over a 
vertical line in the complex plane. 

We recall the elementary formulas 


(5.7) 
and 

(5.8) 


r°° dt 

T{w)a~ w = / e~ at t w —, Reu; > 0, 
Jo t 




k =0 


(~i) A 

k\ 


Let fm,i (t) be the functions defined in (5.4). By using (5.5) and (5.6) we compute the Mellin 
transform of the product of these functions. Observe that the asymptotics for the modified Bessel 
function ( 8 . 8 ) and (8.10) ensure that the Mellin transform f mi fm 2 {s) is well defined for 0 < s < 1. 
We write 


|r(-s)|/F s (m) = fmjm 2 {s) = / f mi (t)fm 2 (t) 


fm 2 (t) 


2iri 


dt 

P+~ s 

/»Cl+ZOO 


f mi (xi)t Xl dx\ 


Cl~too 


dt 

t l+s 


pci+zoo 


fm 2 (t) 


dt 


L -/0 


fl + S — Xl 


~2ni frnAxi) 

6 j Cl—ZOO 

^ rci+ioo _ _ 

= X —- frm (%l)fm 2 (s — xi) dx\. 

Z7TZ Jd—ioo 


dx i 


Note that the integrals that define f mi (x i) and f m2 (s — x\), see (5.5), are absolutely convergent for 
—mi < — Rea;i < 1/2 and —m 2 < Rexi — s < 1/2, respectively, with Rexi = c\. Thus, it is justified 
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that we can write each f mi above as the inverse Mellin transform of / mi , with c, £ (0,1/2). With the 
computation above, and by (5.5), (5.4) and (5.3), we have that 


(5.9) fmjm 2 (s ) = 


F(mi — X\) r(m 2 +a;i — s) 


27 r 2 i 


/*Ci+200 

/ r(l/2 + a;i)r(l/2-Xi + s) ~ — :: // L rT " v '^ ' pn dxi ' 

/ ci —zoo T(toi + aq + 1 ) r(m 2 — aq + s + 1 ) 


Observe that we have an analogous expression just by interchanging the roles of mi and to 2 . 
By (5.2), when mi,m 2 —» oo, we can write 

r(mi-iri) r(m 2 +a,’i —s) 

r(mi + xi + 1 ) r 

where 


r(mi + xi + 1 ) r(m 2 - Xi + s + 1 ) 
1 „ 1 


= / + // + /// + IV, 


I := 




II := 


III := 


TO 


1 


—2a?i+2s+l 


m. 


o 


1 

2a?i+3 J 5 


m 


IV := O 


2 x !+1 


o 


1 


— 2a?i+2s+3 J ' 


m 


m 


i 

2xi+3 


o 


i 


TO 


— 2x i —|—2 s —|—3 


The term I is the leading term. We compute the integral in X\ containing this term. In order to do 
that, we move the contour of integration to infinity (in the negative direction of the Reaq-axis) and 
then we use the residue theorem (the integrand has poles at X\ = —j — 1/2 with residues equal to 

(-i)ViQ: 

/»ci+ioo 


S:=^~ 

2m 


— 1—2tci — 1—2s+2tci 


m 


T(l/2 + aq)T(l/2 — X\ + s) dx\ 


»rW"i:(3) 


2y'+l/2(_l)J 


7=0 


j'- 


T{l+j + s) 


= 777,1 


oo 

'£( 

7=0 


TO 


2\ j+l+s (-1)7 


j! 


T(l+j + s) 


= mr 2 - 2s £ 

j= o 1 ' 0 


where, in the last equality, we used (5.7) with w = 1 + j + s and a = —J. Observe that, by (5.8), we 
have 


S = to 1 


+i t 


t s dt = T (1 + s)mi^~ za (—1 + 1 


to| 


TOT 


-(l+») 


r(l + a) 

(mf + to|) 1+; 


■ J- 


We obtain that the integral coming from (5.9) involving I is 


4 S T(1 + s) 


12 +2s 


It can be proved that for II, III and IV (the error terms) we have better decay in to. □ 

Remark 5.2. Lemma [577] can also be proved for any dimension n > 2. The proofs of these statements 
need rather cumbersome computations, inspired by [5) Subsection 6.4], with suitable modifications, 
and we omit them for the sake of clarity. 

Remark 5.3. In order to prove analogous “asymptotic estimates” for the kernel of the fractional 
discrete integral by using Mellin transforms, in two dimensions, one ends up to an expression of the 
form 


r°° rlf 

J e- 4 t J mi (2t)7 m 2 (2t)^ 


4 s 1 


pci+ioo 


= -ir-.\ r(l/ 2 -aq)T(l /2 + aq-s) 


T(m 1 +x 1 ) T(to 2 — aq + s) 

T(to i — xi + 1 ) r(?n 2 + Xi — s + 1 ) 


dx\. 
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There is a discussion about this identity when m\ = m 2 = 0 in 0 Section 6.2] that can be of interest. 
Thus in a similar way to obtain 

(5.10) 0 < A'_ s (m) = — 2 ' 2 1 2 —f higher order terms, when both m* —> ± 00 , i = 1,2. 

IHI2 S 


6. Illustrations in dimension one 


As we said in the Introduction, we may raise the question whether the Poisson problem (—A ) s u = / 
can be discretized by using the formulas we obtained for the fractional discrete Laplacian and fractional 
discrete integral. One would like to see if the solutions to the discrete Poisson problem converge in 
some sense to the solutions of the continuous Poisson problem. 

Some numerical experiments by using a numerical method were already made in m- Our intention 
is not to make numerical experiments here, but to illustrate the convergence of the discretized problem 
to the continuous one. Furthermore, with the analogous formula in one dimension for the fractional 


discrete integral kernel given in (4.1) we are also able to solve the discrete Poisson problem, so we 


improve in a sense the examples performed in [Ill- 

Formula (1.3) provides an exact expression for the kernel of the fractional discrete Laplacian in 


dimension one. This and the asymptotics in that formula will allow us to compute and draw some 
approximations in one dimension. 

We present some pictures of known examples in dimension one. Some of the examples included in 
this section and in Section [7] were studied, with different numerical methods, by Huang and Oberman 
in mj. 

By using © and O and by taking into account (|5.2[) we can write 


(-A h ) s Uj « F\(j, h, s ) + F 2 (j, h, s) + F 3 (j , h, s), j £ Z, 


where 

( 6 . 1 ) 




E r(H_S) 

|m| <iV 
m^O 


r(|m| + 1 + s) ’ 


( 6 . 2 ) 

and 

for certain N £ N. 


F 2 {3,h,s)~^ s u j E 


1 


|l+2s 


|m| > TV 


F 3(j,h,s):— 2g E u 0 -m I |i_|_ 2 a ’ 

|m|>A 1 1 


On the other hand, formula (4.1) allows us to solve the discrete Poisson problem for a given datum 


/. By using (4.1) and again noticing (5.2) we can write 


(~A h ) s f j =U 1 (j,h,s) + U 2 (j,h,s), j £ Z, 


where 

(6.3) 

and 


Ui(j,h,s) := c_ s h 2s E fj-r 


|m|<IV 


r(|m| + s) 
T(|ro| + 1 — s) 


U 2 (j,h,s) := c- s h 2s E fj-r 


1 


11 — 2 3 ’ 


|m| >7V 


for certain N £ N. 
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6.1. Example 1. The first example to be examined is the function 
(6.4) u(x)=e~ x2 . 

The fractional Laplacian of u at x = 0 can be obtained exactly by using Fourier transform, 

4 s T(l/2 + s) 


(6.5) 


(—A) s m(0 ) = -= / k 2s e~ k2/i dk = 
V 7r Jo 


see m Section 6.1]. Since the function u(x) in (|6.4[ ) has rapid decay, we may ignore the term 
F 3 (j. h, s). In Figure[l]we can see the exact value of (6.5) and the approximation F] (j, h, s)+F 2 (j, h, s ) 
in (16.11) and (16.21) for the fractional discrete Laplacian related to r^it. 


1.0 




• 

• 

0.8 

• 

• 

• 

• 

• 

0.6 

. 

0.4 

• 

• 02 

• 

• 

• 

-2 - -1 

1 * . 2 


• 



-0.2 

L 


Figure 1. The small-dotted line represents Fi (j , h, s ) + F 2 (j , h, s ) for the function 
u(x) in (6.4), with N = 1000, s = 0.25, h = 0.1 and j £ Z, —20 < j < 20, so that the 
horizontal axis is [— 20h, 20h]. The blue point is the value in (|6.5|). 


6.2. Example 2. The second example is the function 
(6.6) u(x) = (l + x 2 )- (1/2 " s) 

with the exact fractional Laplacian 


(6.7) 


(-A ) s u(x) 


4 s r(l/2 + s) 
F(l/2 — s) 


(1 + x 2 ) 


(l/2+s) 


fix), 


see nn formula (41)]. 

The function u(x) decays algebraically, and we are going to ignore the term F 3 (j,h,s) again. 
Actually, it is possible to take into account that this decay is as (/im) - ^ 1-2 ^ to estimate a value for 
F 3 (j, h, s). Some comments and strategies about this are given in |11| . But our goal in this section 
is to illustrate our results, so we will not worry about it. Moreover, we are using N = 1000 so F 3 
is not as important as it would be for small values of N. In Figure [2] the continuous line represents 
f(x) in (6.7) and the dotted line is the approximation F\{j,h, s) + F 2 (j,h, s) in (6.1) and (6.2) for the 
fractional discrete Laplacian of r^rt, with u as in (6.6). 

















FRACTIONAL DISCRETE LAPLACIAN VS DISCRETIZED FRACTIONAL LAPLACIAN 


15 


0 . 20 . 


T).I5 - 


0.10 


0.05 



-4-2 24 


Figure 2. The dotted line represents F\ (j, h, s) + F 2 (j , h, s ) for the restriction of the 
function u(x) in ( 6 . 6 ), with N = 1000, s = 0.4, h = 0.1 and j £ Z, —50 < j < 50, so 
that the horizontal axis is [—50ft, 50ft]. The continuous line is f(x) in (6.7). 


For the discrete Poisson problem, in Figure [3] we plot Ui(j,h,s) in (|6.3[) for r^f, with f{x) as in 


(6.7). We omit the term U 2 {j, ft, s ) because of the decay of the function /( x) (we again remit to [TT] 
for strategies to take into account this term). The continuous line is u(x) in ( 6 . 6 ). 



Figure 3. The dotted line represents U\{j,h,s) for the function (6.7), with N = 
1000, s = 0.4, h = 0.1 and j £ Z, —50 < j < 50, so that the horizontal axis is 
[—50h,50h], The continuous line is u(x) in ( 6 . 6 ). 


6.3. A preliminary lemma. Before continuing with the next examples, we first prove the following 
lemma, see also [I] Lemma 2.4] (note that there are some typos therein). Analogous computations, 
but for the fractional Laplacian instead of the fractional integral, are carried out in [5]. Nevertheless, 
we include our own proof for completeness. Since the lemma is stated for any dimension n > 1, we 
will denote ||®|| 2 := (£” = i \xi\ 2 ) 1/2 , for 2 £l" 


Lemma 6.1. Let 0 < s < 1, 7>0 and x £ K™. Let 

f(x) = (i - ll^lll)X /2 , 

where we use the notation = max{y, 0}. Then 


c. 


tp ( n—2s 


7 ,s,n 2 


C lt s,n\\x\\l S ~ n 2 F 1 


( n— 2s 2—2s . 


)■ 


n +7 

2 



u(x) := (—A) s f(x) 


if IM| 2 < 1, 
if Iklb > 1, 
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where 


and 




- 2a r(^)r(2 + i) 


r( 


2s+7 


i)r(f) 




- 2 s r(^)r(i + i) 


r(^±2 + i)r(s) ’ 
and 2-F\ is the Gaussian or ordinary hyper geometric function. 

Proof. On one hand, it is known that (see [251 P- 171]) 

■77(0 = ^' 7/2r (i +1) neii2 2r ^n +7 )/ 2 (2 7 r|ieii 2 ), 

where J a is the Bessel function of order a. On the other hand (see for example m ch. V.l]), 

u(x) = (-A)- S f(x) = (2n)-^T-\\\ ■ \\f 2s Tf)(x). 


Since / is a radial function, the Fourier transform of / is in fact the so-called Hankel transform 7if, 
and H -1 = TL. Therefore (see [25J P- 155]), 

(6.8) u( x) = r (h + l ) r lleil 2 _2S_i ^ +7 )/2(27r||^|| 2 )J n/2 _ 1 (2^||C|| 2 ||x|| 2 ) 

22s- 1 || a .||| / J o 

We will use the identities (see [25] Ch. 13.4, p. 401]) 


(6.9) 


and 


( 6 . 10 ) 


Jp(at)J v (bt) 

T x 


dt = 


b v r(t 


]^iypf fL+is—X+l a 
2 > 


2 x a v ~ x+1 T(iz + l)r( A+M ~ 1/+1 ) 

p + v — A + 1 v — A — u + 1 h 2 

;^ + i; 

a 


X 2-F 


b 2 \ 

; ^ + 1; , 0 < 5 < a, 

a z / 


Jfj,(at)J„(bt) 


dt = 


a M r( ^~ A+1 ) 


2 A 6 ^- A +iT(/i + i)r( A+t/ ~ M+1 ) 


X 2 F;l 


p -v — A + 1 /i — A — p T 1 


a 2 \ 

0 < a < 6, 


2 2 

which are valid when p + v — A > — 1 and A > —1. Now, by taking p = £±2, v = % — 1, A = 2s + 
a = 2-7T and b = 27r||x|| 2 , we apply (6.9) in the case ||x|| 2 < 1 and (6.10) in the case ||x|| 2 > 1 to 
compute the integral in (6.8). The desired result follows. □ 

6.4. Example 3. The following example involves the function with compact support 


( 6 . 11 ) 


f{x) = {l-x) 


2\ 1—s 


By Lemma 6.1 the solution of (—A ) s u = f is given by 

^^Jjjjg4- s r ( i/2 -,)r( 2 -,) (l _ (1 _ 2a )a; 2 ), if |*| < 1, 

2^1 (l/2 - s, 1 - s; 5/2 - s; ^), if |z] > 1, 


( 6 . 12 ) 


u(x) = 


see also [T] and m (42)]. 

In Figure [4] the dotted line represents the sum Fi(j, h, s ) + F 2 (j, h, s) in (6.1) and (6.2) related to 
the function u(x) in (|6.12|). 
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Figure 4. The dotted line represents Fi(j, ft, s) + F 2 (j, ft, s ) for the function (6.12 ), 
with N = 1000, s = 0.25, ft = 0.1 and j £ Z, —20 < j < 20, so that the horizontal 
axis is [— 20ft, 20h\. The continuous line is f(x) in (6.11). 


The continuous line is f(x) in (6.11). For solving the discrete Poisson problem in this case, in 
Figure[5]the dotted line represents Ui(j,h,s) in (6.3) for the function f(x) in (6.11). Since /( x) has 
compact support, we can choose N such that the term U 2 (j, ft, s) is zero. The continuous line is f(x) 
in ( |6.12[ ). 



Figure 5. The dotted line represents Ui(j, ft, s ) for the function (6.11), with N = 20, 
s = 0.25, h = 0.1 and j £ Z, —20 < j < 20, so that the horizontal axis is [—20 h, 20ft.]. 
The continuous line is u(x) in (6.12). 


6.5. Example 4. We present another function with compact support: 
(6.13) /Or) = (1 - z 2 ) 2 + - s . 


By Lemma 6.1 the solution of (—A ) s u = f is given by 

V s r(1/2 ;?£ (3 ~ s) (l - (2 - 4s)x 2 + (!-§«+ | S 2 K), 


(6.14) 


i(x) = 


if |x| < 1, 


k 4 ~ s r r('l)r(7/2-7) S) kl 28 " 1 2^1 (l/2 - 8 , 1 - s ; 7/2 - a; ^) , if |*| 


> 1 , 


see also [T] and m formula (43)]. 

In Figure[6]the dotted line is F] (j. ft, .s) + F 2 (h,j,s) in (6.1) and (6.2) for u(x) as in (6.14). The 
continuous line is the exact F(x) in (6.13). 
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Figure 6 . The dotted line represents Fi(j, h, s) + F 2 (j, h, s ) for the function (6.14), 
with N = 1000, s = 0.25, h = 0.1 and j £ Z, —20 < j < 20, so that the horizontal 


axis is [— 20h, 20h\. The continuous line is f(x) in (6.13). 


For the discrete Poisson problem, in Figure[7]the dotted line is Ui(j, h, s) in (6.3) for (6.13). Again, 
since f(x) has compact support, we can choose N large such that U 2 (j,h,x) is zero. The continuous 
line is u(x) in (6.14). 



Figure 7. The dotted line represents U-\ (j. h, s ) for the function f{x) in (6.131 , with 
N = 20, s = 0.25, h = 0.1 and j £ Z, —20 < j < 20, so that the horizontal axis is 


[—20/i,20/z] . The continuous line is u(x) in (6.14). 


7. Illustrations in dimension two 

In this section we plot approximations of some known examples in dimension two. 

Unlike the one dimensional case, we do not have explicit expressions for the two dimensional kernels 
of the fractional discrete Laplacian and discrete integral. Therefore, computation and evaluation 
of the kernels is complicated, so we are going to use a less precise fast method. We only have 


at our disposal the asymptotic estimates for the kernels K s (m) and /\_ s (m) of Lemma 5.1 and 


Theorem 4.1 Nevertheless, in view of the one dimensional case and Lemma 8. 2[ we conjecture that 


we can approximate the kernel of the fractional discrete Laplacian by the main term in (5.1), and the 


kernel of the fractional discrete integral by the main term in (5.10). In this way, for the fractional 
discrete Laplacian we take 




mGZ \{0} 


VLj 'LL n — u 

"-il 2 +2s 


+ F —: F(j, h 1 s) + E, j £ Z 2 


(7.1) 
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for certain N £ N, where E involves the sum with the error terms. 


On the other hand, we use formula (5.101 as an approximation that allows us to solve the discrete 

Jj—m 


Poisson problem for a given datum /. We write 


(7.2) 


= (-A h )-°f j = C- 2 , s h 2 °(f j K_ s ( 0 )+ E 


E 


mGZ 2 \{ 0} 
I | < AT 


U(j,h,s) + E, j £ Z 2 , 


for certain N £ N, where E involves the sum with the error terms. Moreover, A'_ s (0) in (7.2) can be 
explicitly evaluated. Indeed, 

K - M = 4 r G(0 - ‘>4 = r hf e " 4 ‘ (/ ” (a,)2 4 = 4 " ^ 44-4 I; 

where aiy is the generalized hypergeometric function. 


7.1. Example 1. The first example we deal with is a function with compact support: 

(7.3) /(*) = ( 1-IMI2&-, 


By Lemma 6.1 the solution of (—A ) s u = / is given by 


(7.4) 


(4-T(l- S )r(2- S )(l- (1-5)M1), 

U '- X ) | 4—s r (!—g) r (2—s) 11 ||2s—2 p A 1 o. Q J- 1 

r(«)r( 3 -s) \\ x \\2 2 T 111 s, i s, o s, || x ||; 


if IMI 2 < 1, 
if \\x\\ 2 > 1. 


In Figure [8| we plot the exact function f(x) in (7.3) on the left, the approximation F(j, h. s) in (7.1) 
related to (7.4) in the center and the error (difference) on the right. 



Figure 8. The surface on the left is f(x) in (7.3), with x £ [—2,2] x [—2,2]. In the 
center we find the approximation F(j, h, s) related to the function u(x) in (7.41, with 
N = 500, s = 0.25, h = 0.1 and j = {ji,j 2 ) £ Z 2 , —20 < ji < 20. The mesh is 
[—20/i,20/i] x [—20/i, 20/i]. The error in each point of the considered mesh is plotted 
on the right. 


In Figure [9] w e plot the exact function u(x) of (7.41 on the left, and the picture arising from 
U(j, h, s ) in ( |7.2[ ) related to (7.3) on the center. This situation corresponds with the discrete Poisson 
problem. The error is drawn on the right. Observe that, since / is a function with compact support, 
there is only a finite number of nonzero terms in the computation of U(j, h, s). 
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Figure 9. The surface on the left is u(x ) in (7.41 , with x £ [—2, 2 
center we find the approximation U ( j , h, s ) related to /(x 
s = 0.25, h = 0.1 and j = (ji, j 2 ) 


x [—2, 2]. In the 
in (17.31) , with N = 40, 


£ Z 2 , —20 < ji < 20. The mesh is [—20/i, 20/i] x 
—20/i, 20/i]. The error in each point of the considered mesh is plotted on the right. 


7.2. Example 2. The second example concerns the compactly supported function 


(7.5) 


f(x) = (i-\\xr 2 y + -°, XG 


By Lemma 6.1 the solution to (—A) s u = f is given by 
(7.6) u(x ) = 


4 -s 2 -1 r(l - s)T(3 - s)(l - (2 - 2s)||m||l + (!-§«+ ^ 2 )||x||t), if ||x|| 2 < 1, 


4_s TWpir ll^ll' 8-2 2*1 (l - s, 1 - s; 3 - s; 


if IIccII2 > 1. 


In Figure [To] we see the picture of f{x) of (7.5) on the left, the approximation F(j, h, s ) in related 
to (|7.6|) in the center, and the error on the right. 



Figure 10. The surface on the left is f(x) in ( |7.5| ), with x £ [—2,2] x [—2,2]. 
In the center we find the approximation F(j,h,s ) related to u(x) in ( |7.6| ), with 
N = 500, s = 0.25, h = 0.1 and j = (j‘1,.7‘2) £ Z 2 , —20 < ji < 20. The mesh is 
[—20/i,20/i] x [—20/i, 20/i]. The error in each point of the considered mesh is plotted 
on the right. 


Next we illustrate the problem of solving the discrete Poisson problem in this case. In Figure 0 
the picture on the left is u(x) given in (7.6), the approximation U(j,h,s) in (7.2) related to f{x) in 
(7.5) is in the center (observe that there are finite summands, so we can choose N large so that the 
term E is zero), and the error is on the right. 
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Figure 11. The surface on the left is u(x) in (7.6), with x £ [—2,2] x [—2,2]. 
In the center we find the approximation U(j,h,s) related to f(x) in (7.5), with 
N = 40, s = 0.25, h = 0.1 and j = (ji, ji) £ Z 2 , —20 < j,, < 20. The mesh is 
[— 20h,20h] x [—20h,,20h]. The error in each point of the considered mesh is plotted 
on the right. 


7.3. Example 3. The last function we take into account is the following. 
0 < a < 2 — 2 s. Consider the function 

P 2 


Let 0 < s < 1 and 


(7.7) 

It is known that 

(7.8) 


u{x) = 


X 


i 2 ~ q > 


X £ . 


(— A) s u(x) = 2 


2 s T(<r/2 + s)T(l — a/2) 
T(1 — a/2 — s)T{a/2) 


| 2 -“- 2s 


=: /(:r). 


Observe that f(x) in (7.81 has a singularity at the origin. However, we are going to check the behavior 


of the numerical method by computing the corresponding approximations F(j,h,s) in (7.1) for the 


fractional discrete Laplacian and U(j, h , s) in (7.2) for the fractional discrete integral. In order to avoid 


the singularity in the origin, we move the mesh a distance h/2 in the direction of each coordinate 
axis, so we consider a mesh [—(J + 1/2 )h, (J + 1/2 )h\ x [— (J + 1/2 )h, (J + 1/2 )h}. Then, we see in 
Figures 12 and 13 respectively that these approximations are close to the exact functions f(x) in (7.8) 
and u(x) in (7.7), and the error of the approximation is large only near the origin, which is expected 
because both u and / tend to infinity close to zero. 



Figure 12. The surface on the left if f{x) in (7.8) with a = 0.5 and s = 0.3, for 
x £ [—2,2] x [—2,2]. In the center we find the approximation F(j,h,s) related to 
u(x) in ( |7.7| ), with N = 500, h = 0.1 and j = (j'1,/2) £ Z 2 , —20 < ji < 20. The mesh 
is [—(20 + l/2)/i, (20 + l/2)/i] x [— (20 + l/2)/i, (20 + l/2)h]. The error in each point 
of the considered mesh is plotted on the right. 
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Figure 13. The surface on the left is u(x) in ( |7.7[ ) with a = 0.5, for x £ [—2,2] x 
[—2,2]. In the center we find the approximation U(j,h,s ) related to f(x) in ( |7.8| ), 
with N = 500, s = 0.3, h = 0.1 and j = (jujz) £ Z 2 , —20 < j t < 20. The mesh is 
[—(20 + 1/2 )h, (20 + 1/2 )h] x [—(20 + 1/2 )h, (20 + 1/2 )h\. The error in each point of 
the considered mesh is plotted on the right. 

8. Technical lemmas and properties of Bessel functions 

8.1. Some technical lemmas. Lemmas in this subsection are needed in the proof of Theorem o 
They are also useful to get estimates for the kernels of the fractional discrete Laplacian (2.2) and for 


the fractional integral kernel in Theorem 4.1 


The following lemma is elementary. We include it for the sake of completeness. 
Lemma 8 . 1 . Let A > 0. Let a,b be real numbers such that 0 < a < b < oo. Then 


min{A,1} < 


6 A - a x 

6 A_1 (6 — a) 


< max{A, 1}. 


b x — a x 


1 - (a/by 


Proof. Let us first suppose that A > 1. Then 

0 <a< 6 <oo=> 0 < a /6 < 1 => 0 < (a/b) x < a/b < 1 — 

6 A_1 (6 - a) 1 - a/b 

On the other hand, by applying the mean value theorem to the function x x , we get 


> 1 . 


6 a — i 


1 ~(a/by 


= Xx 


A—1 


< A, 


6 A_1 (6 — a) l-a/b 
for certain x £ (a/ 6 ,1). In the case 0 < A < 1, the proof is analogous. 
Lemma 8.2. Let 0 < s < 1, t £ R, and m £ Z, m ^ 0. Then, 


□ 


( 8 . 1 ) 

Also, if 0 < s <1/2, we have 

( 8 . 2 ) 


r(H - s) 


l 


r(|m| + 1 + s) 
r(|m| + s) 


I l+2s 


c. 


T(|m| + 1 - s) 


m 


< 


m | 2 -2s 


Proof. Let us begin with the proof of (8.1). Without loss of generality, take m £ N. We write 


T(m — s) 

1 

< 

T(m — s) 

1 

+ 

1 

1 

T(m + 1 + s) 

m !+2s 

T(m + 1 + s) 

(m — s) 1+2s 

(m — s) 1+2s 

m 1+2s 


The second term can be easily estimated, just by applying Lemma 
6 = 1 , namely 

m— 7 J 


5.1 


with A = 1 + 2s, a = — and 


1 1 

~ 1 

( 1 M 

~ c ' 

(m — s) 1+2s m 1+2s 

(to — s) 2s 

\m — s m ) 

m 2+2s 


m—s 
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where the symbol ~ means that constants depend only on s. Now we study the first term. For k £ N, 
we have (see for instance EZl Section 7]) 

r (k - s) 1 


(8.3) 

With this, 

r(l + 2 s) 


T(k + n + s) r(n + 2s) 
r(m — s) 1 


/*oo 

/ e - {k - s)v (l - e~ v ) n+2s - 1 dv. 

Jo 


r(m + n + s) (m — s) 1+2s 


< 


0 —(m—s)v L,2s 


v 2s - (1 - e~ v y s \dv 


e -(m-s)v v 2s 


e -(m-s)v v 2s 


1 - 


1 - 


(^rT 


dv 


1 — e 


dv 


< 


r e -(--) v^o»~ r (i +2,) 1 

L 2 


where we applied Lemma 8.1 and in the last inequality we used that ^ > u ^ 1 + e for u £ (0, oo). 


The proof of (8.2) is analogous, with the restriction 0 < s < 1/2 coming from Lemma 8.1 


□ 


8.2. Properties of Bessel functions Ik. We collect in this subsection some properties of modified 
Bessel functions. Let Ik be the modified Bessel function of the first kind and order k £ Z, defined as 

2 m-\-k 

(8.4) 4 (f) = r k j k (it) = Y, 


l 


m =0 


to! T(to + k + 1 ) \2 


Since k is an integer and l/r(n) is taken to be equal zero if n — 0, — 1,— 2,..., the function Ik is 
defined in the whole real line. It is verified that 

(8.5) J_ fe (t) = 4(f), 

for each k £ Z. Besides, from (8.4) it is clear that 4(0) = 1 and 4(0) = 0 for k ^ 0. Also, 

( 8 . 6 ) I k {t) > 0 
for every k £ Z and t > 0, and 

J>- 2t 4(2f) = 1. 
fce z 


E e- 2nt Y[i ki m = i. 


From this it can be verified that 

(8.7) 

/cGZ 71 i =1 

On the other hand, there exist constants C, c > 0 such that 

( 8 . 8 ) ct k < h{t) < Ct k , as t —f 0 + . 
In fact, 

k 


(8.9) 


hit) 


1 


for a fixed k ^ —1, —2, —3,... and t —>• 0 + , 


2J T(k + 1)' 
see m- It is well known (see [J3]) that 
(8.10) l k {t) = CeH~ 1/2 + R k (t), 

where 

|4fc(f)| < Cfce^ -3 / 2 , as t —f oo. 
We also have (see [16]) that, as v —> oo, 


( 8 . 11 ) 


hiz) 


1 / \ 1 
1 ez\ 


2 v) 2 u v\' 























24 


O. CIAURRI, L. RONCAL, P. R. STINGA, J. L. TORREA, AND J. L. VARONA 


For the following formula see [TT] p. 305]. For Rec > 0, — Re v < Re a < 1/2, 


( 8 . 12 ) 



e~ ct I v {ct)t a ~ 1 dt = 


(2c)~ a r(l/2 — a)r(a + u) 
y/n Y{y + 1 — a) 
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